Intrinsic spin Hall effect in graphene: Numerical calculations in a multi-orbital model 
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We study the spin Hall effect (SHE) in graphene using a realistic multi-orbital tight-binding 
model that includes the atomic spin-orbit interaction. The SHE is found to be induced by the 
spin-dependent Aharonov-Bohm phase. In the metallic case, the calculated values for the spin Hall 
conductivity (SHC) are much smaller than the quantized Hall conductivity for realistic parameter 
values of metallic graphene. In the insulating case, quantization of the SHC is violated due to the 
multi-orbital effect. The present study suggests that the SHE in a honeycomb lattice is enhanced 
by chemical doping, such as the substitution of carbon atoms with boron atoms. 

PACS numbers: 73.43.-f, 72.25. Hg, 73.61.Wp, 85.75.-d 



The recent discovery of graphene has stimulated much 
interest in electron transport in graphene due to its un- 
conventional electronic structure. [l| The crystal struc- 
ture of graphene is a two-dimensional honeycomb lat- 
tice of carbon atoms. In the absence of the spin-orbit 
interaction (SOI), the band structure of graphene is de- 
scribed in terms of massless Dirac fermions at the edges 
(K and K' points) of the Brillouin zone. Due to a non- 
zero Berry's phase (a geometric quantum phase) Q, a 
striking "half-integer" quantum Hall effect is realized 
in graphene [1, Q|. The nature of the quantum Hall 
effect near the massless Dirac point has been studied 
theoretically. [!, E3, 0, @, @, E3] Concerning transport phe- 
nomena in graphene, the emergence of a sizable spin Hall 
effect (SHE), which is a phenomenon that a spin current 
flows perpendicular to an a ppl ied electric field, had been 
predicted theoretically. [Ill, [l2| In particular, Kane and 
Mele[ll[ have demonstrated that the SOI generates an 
energy gap at the Dirac point, and the "quantum SHE" 
is expected to appear in insulating graphene. 

Recently, the theory of the intrinsic SHE caused by 
the Berry phase in zincblend semiconductors [l3} or by 
a uniform SOI such as the Rashba SOI in a two- 



dimensional electron gas[14( has been attracting much 
interest. Quite recently, the SHE has also been ob- 
served in the simple metal Al[15l] and the transition metal 
Pt.[l(| The spin Hall conductivity (SHC) reported for 
Pt is 240(ft/e)(£lcm) _1 at room temperature, which is 
much larger than that of semiconductors. The huge 
SHC in Pt has been explained theoretically by Kontani 
et al. in terms of the atomic SOI and a realistic tight- 
binding (TB) model that includes inter-orbital hopping 
integrals. [l?! EH, [3] They further showed that the so- 
called current vertex corrections (CVC) have little effect 
on the SHC in transition metal, in contrast to that in 
semiconductors with the Rashba SOI.jiol] 

The successful reproduction the experimental values of 
the SHC in transition metals by the TB model with the 
atomic SOI encouraged us to apply this model to other 
systems. Since the electronic structure of graphene is well 



reproduced by the realistic TB model, we use this model 
to calculate the SHCs in both undoped and doped (metal- 
lic) graphene, and compare the results with those calcu- 
lated using the Kane-Mele model, [ll) which contains only 
a p z orbital and an effective SOI to produce an energy 
gap at the Fermi energy. 

In the present paper, we adopt a realistic multi-orbital 
(s, p x , p y and p z orbitals) TB model with the atomic 
SOI on a honeycomb lattice, and use the Kubo-Streda 
formula 2l| to calculate the SHC. It should be noted that 
the p x and p y orbitals contribute to the SHC because they 
are mixed with the p z orbital via the atomic SOI. By 
calculating the SHC as a function of the Fermi energy, 
we find that the SHC becomes large in the energy re- 
gion where the p x and p y orbitals are dominant, far away 
from the Dirac point. The effect of the CVC, which is 
calculated in the self-consistent Born approximation, is 
appreciable and causes the SHC to double compared to 
that without the CVC, while the qualitative behavior of 
the SHC remains unchanged. In the case of insulating 
graphene, SHC is not quantized because s z is not con- 
served due to the SOI between the p x {py) and p z orbitals. 

The Hamiltonian for electrons on the honeycomb lat- 
tice, which is decomposed into A and B sub-lattices, is 
given as H = H +Hso where Hq and Hso are the kinetic 
and SOI terms, respectively. Hq is given as 
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where t"? denotes the nearest-neighbor hopping integral 
between the a orbital at site i and the (3 orbital at site 
j, and a denotes the spin. The Slater-Koster parameters 
in the TB model are taken as ssa — —0.86, spa = 1.14, 
ppa = 1.0 and pp-n = —0.5 in units of ppg(5.7eV).[22| 
Hereafter, we take ppa and the lattice constant a(2.55A) 
as the units for energy and length, respectively. We also 
put H = 1 . 

Hso represents the atomic SOI given as 
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where A is the interaction constant, and I and s are the 
orbital and spin angular momentum, respectively, -ffso 
has off-diagonal elements such as (p x ± |i?sobj/±} = 
=rA/2, (p x ± \Hso\PzT) = ±A/2 and (p y ± l-ffsol^T) = 
—iX/2, where \la) represents an electron with orbital I 
and spin a = ±1. 

According to Streda, [2l( the intrinsic SHC at T = is 



given as crjL = a z x \ + erf 1 , where 
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Here, G R ^ is the retarded (advanced) Green's function 

given as G R(A \k,uj) = fu> — H + (— )iT\ , where f is 
the matrix form of the imaginary part of the self-energy 
(damping rate) due to scattering by local impurities. The 
Green's functions are represented by 16 x 16 matrices in 
the momentum representation. 

The matrix forms of the charge and s z -spin current 
operators for the /u-direction (/i = x, y) are given by 

= ~ e W; and ^1 = respectively. In 

the SHC, a xy and er** 1 represent the "Fermi surface term" 
and the "Fermi sea term" , respectively. In this paper, we 
take N = 900 x 900 fc-point meshes in the numerical 
calculations. 

We first neglect the CVC and assume that the damp- 
ing matrix is diagonal and is independent of orbital 
= 7<5 Q/ 3. Figures Q] (a) and (b) give the calculated 
results of the SHC as a function of the position of the 
Fermi energy Ep and the band structure of graphene for 
A = 0.1 and 7 = 0.01 (in units of ppcr), respectively. In 
Fig. [2(a), the solid and dotted lines denote a xy and a x l y , 
respectively. We observe that the SHC a xy is large in en- 
ergy regions where the bands are nearly degenerate, while 
a x l y exhibits peak structures in these energy regions. In 
the metallic region, the magnitude of a x l y is much larger 
than that of a^ y l . 

The insets of Figs. QJa) and (b) show the SHC and 
the band structure near the Dirac point of graphene, i.e., 
Ep = 0. The Dirac point at K is split by the atomic 
SOI, and forms an energy gap with Ad = 0.0013 when 
A = 0.1. 

In order to study the SHE in graphene, we estimate the 
value of 7 to be about 0.005 by comparing the longitudi- 
nal conductivity a xx calculated in the present model with 
the experimental resistivity observed for doped graphene 
p = l/cr xx = 10017. [2^1 In this case, the number of 
electrons n per atom of doped graphene is n = 3.998, 
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FIG. 1: (a)SHC af v + crfy (solid line) and afy (dotted line), 
and (b) band structure through k points K= (0,47r/3) , T = 
(0,0) and M= (7r/v3,7r) against energy for atomic SOI A = 
0.1 and damping 7 = 0.01 (in units of ppa). The insets show 
the SHC and band structure near the Dirac point. 



which corresponds to Ep = —0.05 in the present model. 
The value of A of the atomic SOI at the Dirac point in 
graphene is estimated to be A = 0.001, 24, 2|| which cor- 
responds to an energy gap Ad — 2 x 10~ 7 (~ 10 _6 eV). 
The SHC calculated by using A = 0.001 and 7 = 0.005, 
and by neglecting the energy dependence of 7 is a xy ~ 
3 x 10~ 5 [e/27ra] at the Dirac point. This value is much 
smaller than the quantized SHC e/27ra, since 7 ^> Ad 
and the graphene is metallic for these parameter values. 

The SHC may increase in doped graphene; for ex- 
ample, the value of a xy increases to ~ 0.1[e/27ra] at 
Ep = —0.76, which can be realized by substituting 
50% carbon atoms with boron atoms. Although this 
situation may be virtual, similar electronic state may 
be realized by graphite intercalation. Since l[e/27ra] ~ 
1500(fi./e)(ricm)~ 1 in the present case (a = 2.55A), the 
SHC - 150(ft/e)(ficm)- 1 at E F = -0.76 is the same 
order as that of Pt. 



Now, we give an intuitive explanation of why the SHE 
is induced by the atomic SOI. The calculated A depen- 
dence of the SHC is shown in Fig. [2j We see that 
a xy oc A and a xy oc A 2 for n = 3.2 (Ep = —1) and n = 4 
(Ep = 0), respectively. The results may be interpreted 
as follows. Since the p x and p y orbitals are dominant at 
Ep ~ —1, an anti-clockwise motion of an up-spin elec- 
tron on p y orbitals in a honeycomb lattice may be given 
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FIG. 2: a% y against A with 7 = 0.01 (in units of ppa) for 
n = 3.2 (E F = -1) (solid line) and for n = 4 (J5 F = 0) 
(dotted line), respectively. 



FIG. 3: Effective AB phase in honeycomb lattice derived by a 
motion of an up-spin electron (a) on mainly p x and p y orbitals, 
and (b) on mainly p z orbital. 



by Fig. [3]Ja). The arrows in the figure represent an inter- 
orbital transition induced by the atomic SOL In the pro- 
cess p y — > p x (Px — ► Py), the SOI works once, and yields a 
factor (— )iX/2, which is the first order of A. The corre- 



sponding motion of an electron yields a factor i = e 



2-Ki/i 

: 

which can be interpreted as the Aharonov-Bohm (AB) 
phase factor e 27 ™^/* (<^> = hc/\e\), where 4> is effective 
magnetic flux (f> — § A - dr = 0o/4 through a honeycomb 
lattice. Since the sign of the effective magnetic flux is 
opposite for down-spin electrons, up- and down-spin elec- 
trons move in opposite directions under an electric field. 
By this process, <r xy oc A is realized in the region where 
the p x and p y orbitals are dominant, as indicated by the 
solid line in Fig. [2j 

At Ep ~ the p z orbital is dominant, and an impor- 
tant process of an up-spin electron on p z orbitals may 
be given by Fig. [31(b) . Solid arrows and a dotted arrow 
represent the inter-orbital transitions induced by the SOI 
and by hopping, respectively. The SOI should operate at 
least twice while the electron moves around the honey- 
comb structure, since the spin is flipped in the transitions 
Pz — > Py and p x — > p z via the SOI in this case. The cor- 
responding motion yields a factor of i = e 27 ™/ 4 , which 
corresponds to the AB phase and an effective magnetic 
flux. By this process, a xy oc A 2 is realized in the region 
where the p z orbital is dominant at the Fermi level, as 
indicated by the dotted line in Fig. [2J. Since A <C 1, the 
SHC near Ep — is about 10 times smaller than that 
near Ep = — 1. 

Up until now, the CVC have not been taken into ac- 
count. However, the CVC may play an important role 
for the SHC in graphene, as discussed in Ref. [T3| for the 
Kane-Mele model. In order to study the role of the CVC 
in the present four-orbital model, we employ the self- 
consistent Born approximation, where orbital-dependent 
damping Y afs = j a S a( 3 is given by 
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where ni and I are the density of impurities and impu- 



rity potentials, respectively. The total current (J^ = 

+ AJ^) with the CVC (AJ°) is given by the Bethe- 
Salpeter equation: 



jc(fc,c) = i c (fc)- 
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which is solved self-consistently. Here, we put n-J 2 = 
0.02 , for which the experimental value of resistivity [23| 
p = 1/(T XX = 100S1 at Ep — —0.05 is realized in the 
present calculation. The CVC part for the Fermi surface 
term is obtained from 
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and the total SHC with the CVC is given by a zI y = a zl y + 

The calculated results of a zl y and a zl y near the Dirac 
point are shown in Fig. [4] for A = 0.1, which makes 
Ad = 0.0013 at the Dirac point. The numerical results 
for \E\ < 0.004 are omitted because of poor convergence. 
We note that 5** = in the insulating system (Ep = 0) . 
We see that a* y with the CVC is almost double a x l y with- 
out the CVC. Although the behavior of a z xy is consistent 
with the results of Sinitsyn et al. [12J , a xy is considerably 
smaller than that obtained by Sinitsyn et al. We consider 
that the disagreement comes from the difference between 
two models. The reason why the CVC remains for the 
atomic SOI in a honeycomb lattice is that the Hamilto- 
nian breaks inversion symmetry i? Q ^(fc) ^ H a p(—k). If 
H a p(k) = H a j3(—k) is satisfied, the CVC vanishes iden- 
tically. We note that the CVC is absent for the Fermi sea 
term in the Born approximation in the present model. 

In the Kane-Mele model, the SHE o% y — e/2-K is quan- 
tized at the Dirac point of graphene. [ll| In order to com- 
pare with their results, we calculate a xy in the energy gap 
by taking the limit 7^0, which corresponds to insulat- 
ing graphene. We obtain a x l y = and a zI y ~ \.2e/2ir at 
Ep = for A = 0.4. The reason for the violation of the 
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FIG. 4: Total af y (solid line) with the CVC and af y (dotted 
line) without the CVC for A = 0.1 and A D = 0.0013. 



quantization is that s 7 . is not conserved due to the SOI 
in the present model. 2(| 27 1 

So far we have neglected the effect of lattice deforma- 
tion (curvature effect), which was pointed out to be im- 

" [H, [2^ The curva- 



portant for e.g. carbon nanotube. [2E 
ture effect induces a hopping (At) between p z orbital and 
p x (Py) orbital. At is estimated as ~ 0.01 in graphene [30j] 
and ~ 0.1 in nanotube [3ll] in units of ppa. The lowest or- 
der of the curvature effect on SHC a xy curv is estimated 
as (At) 2 A, where a denominator of order leV is omit- 
ted. Since the SHC a* int without the curvature effect 
is the order of A 2 , the ratio a xy curv /& xy int ~ 0.1 in the 
graphene, and ~ 10 in the nanotube by adopting a re- 
alistic value of A = 0.001 for the SOI of carbon atoms. 
We expect that the curvature effect for the SHC may be 
small in graphene but dominant in nanotube. 

In summary, we have studied the SHC in a two- 
dimensional honeycomb lattice using a realistic band 
structure consisting of s, p x , p y and p z orbitals with the 
atomic SOI. The estimated SHC for a realistic value of 
the atomic spin-orbit coupling A and constant damping 7 
for the metallic graphene at the Dirac point is consider- 
ably small. We predict that the SHC will be large when 

the Fermi level is shifted to E F = -0.76( 4.3eV), 

where the p x and p y orbitals are dominant at the Fermi 
level. In the self-consistent Born approximation, the SHC 
a xy with the CVC is almost double the SHC without the 
CVC. The obtained value of a x l y is considerably smaller 



than that of the Kane-Mele model. 11, 12] In the case of 



insulating graphene, the obtained SHC a xy is not quan- 
tized because s z is not conserved due to the SOI. 
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